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Abstract 


By building on the work in Kuzmak & Tkachuk, “Preparation of 
quantum states of two spin-i particles in the form of the Schmidt 
decomposition”, Physics Letters A, 378, ppl469-1474, which outlined 
the control of the degree of entanglement within this system, it is 
proven that any SU( 4) state manipulation operator can be realised 
for this system using a sequence of pulsed magnetic fields in either 
two linearly independent directions if the gyromagnetic ratios are un¬ 
equal or three directions for equal gyromagnetic ratios. To achieve 
this goal, an elementary Lie theoretic proof of the fact that the group 
of transformations generated by finite products of exponentials of a 
set of Lie algebra vectors is equal to the Lie group generated by the 
smallest Lie algebra containing those vectors is rewritten into an ex¬ 
plicit algorithm. A numerical example as well as the proof of the 
algorithm’s effectiveness is given. 
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1 Introduction: Physical System Overview 

The physical system considered here is the four dimensional quantum state 
space of two, two dimensional qubits. Furthermore, in this paper, the two 
qubits have an unremovable coupling between them, defined by a term pro¬ 
portional to a x ® a x + a y <g) o y + o z ® a z as an interaction term in the 
two-particle system’s Hamiltonian. For concreteness I shall consider a lone 
donor 31 P atom (the naturally-occuring isotope) in a lattice of silicon as in 
[5]. This donor atom’s nucleus has a spin | and at erogenic temperatures the 
lone 31 P in the Si lattice can trap an electron. In this cryogenic, electron- 
trapped configuration, the 31 P donor-electron pair is a coupled-spin, two 
qubit system. With the quantum state space defined by T-L n ® R. ei where 
R n , He are, respectively, the separate two-dimensional spin state spaces for 
the free donor and electron, the Hamiltonian, Schrodinger Equation and full 
coupled spin ket for the coupled system (with the ground state set to energy 
0) in the presence of a magnetic field are: 
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H = Y, {Bj ('YnCTj <g> id - 7e id <g) aj) — § Uj <g> 07) 

j£{x,y,z} 

d tij) = e l&t Tp] i/j d = ( z t Z2 z 3 ^4) 

|V>) '= *1 IV’+i) ® l^+i) + ^2 |^+i) ® IV’-l) + ^3 IV’-l) ® l^+l) + -4 IV'-l) ® l^ll) 

( 1 ) 

where a 3 are the 2x2 Pauli spin matrices, Bj are the Cartesian compo¬ 
nents of the magnetic induction and the subscripts ±1 stand for the spin 
up/spin down kets in the separate nucleus ( n ) and electron (e) spin states. 

7 n , 7 e are the nucleus’s and electron’s gyromagnetic ratio, respectively; these 
are set by the fundamental physics of these particles and cannot be con¬ 
trolled by the experimenter. Their rough values in SI units, as given in [5j, 
are 7 n = 17.23 MHz T ^ 1 and 7 e = 27.97 GHz T _1 . Likewise, the value of 
k ~ 58.765MHz is set by the geometry of the overlap between the electon’s 
orbital and the nucleus: it is always present and cannot be experimentally 
changed. Lastly, there is currently no technology through which an experi¬ 
menter could impart a different magnetic held separately to the nucleus and 
the electron; indeed, significant magnetic held variations over the angstrom 
scales of the system imply lOkeV electromagnetic radiation, which would def¬ 
initely disturb the system in ways not described by the simple model of ([TJ). 

If we make the following dehnitions of the quaternion units and members of 
su( 4 ): 


1 = 





0 

—i 


x 0 = 7-77 ( In i ® 1 + 7 e 1 ® i); Yq = ( 7n j ® 1 - 7e 1 ® j) 

Z 0 = 777 - (- 7 nk<g> 1 + 7 e l <g> k); K = f (i< 2 ) i+ j <g>j + k ( 8 ) k) 

( 2 ) 
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then the physical constraints mean that only time evolution operators of the 
form e x where X is a Lie algebra member of the form 
X = (^B x X 0 + B y Yo + B z Zq ± k K^j r with the SI unit values cited above 
are directly available to the experimenter to control the coupled system’s 
state with. Here r is the magnetic held pulse time. The Bj can be any 
real values within an interval set by the magnetic held pulsing apparatus: 
today roughly — lOrnT < Bj < +10mT is reasonable. The magnitude of the 
co-efficient of K r is always the same and cannot be varied, being roughly 
58.765MHz. However, its sign can be positive or negative, because the one pa¬ 
rameter group {e sA [ s G M} is compact and periodic, i.e isomorphic to U( 1). 
The eigenvalues of K are | i and —\i (the latter a triple), so that e Kt is peri¬ 
odic in t with period 4 7T. Thus, the period of e KAt is Att/k fa 214ns. There¬ 
fore a matrix of the form e~ K x to is of the form exp — to) ^ • e Kt is in¬ 

teresting from the experimenter’s standpoint because it periodically entangles 
then unentangles any given input state, and thus, as shown in [5], a sequence 
of operations of the form e Xr where X = ^ B x X 0 + B y Y 0 + B z Z 0 ± k K^j t 
can be used to prepare any quantum state in B n <g) % e given the input of 
one particular state that is readily prepared by the recipe of [8] . As another 
example, the degree of entanglement is defined in |lj as the von Neumann 
entropy of the renormalised projection of the whole state vector in T-L n ® Be 
onto either Bn, Be, which projection, considered in its subspace alone, is 
equivalent to a mixed quantum state. The degree of entanglement can be 
shown to be a monotonic function of = 2 \zi z 4 — z 2 z 3 \ G [0, 1], both the 
entropy and vary between 0 (for a product state) and 1 (for a maximally 
entangled state) and, for example, e Kl<t { 1 0 0 0 ) T swings periodically 
between a product state and a maximally entangled one. 

In §|3] I shall use Lie theoretic ideas to find a sequence of operations of the 
form e Xr to be equivalent to any unitary operator in S77(4), even though the 
subspace of su(4) spanned by vectors of the form 
X = (B x Xq + By Yq + B z Zq ± k K J r is only four dimensional and there- 
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fore very “small” and “uncomplicated” compared to the full Lie algebra su(4). 
S77(4) is a simple Lie group, so one cannot therefore break its description 
down into a semidirect product of smaller Lie groups, which decomposition 
might simply the goal of realising any any unitary operator in SU(4). How¬ 
ever, one particular, non-normal subgroup of S77(4) which is of interest is 
the subgroup of exponentials of the following Lie algebra. 

Lemma 1. The Lie group exp(D) where: 


0 = span R ({z 1 <g) 1, i ® 1, j ® 1, k <E> 1, 1 <£> i, 1 <8> j, 1 <8> k}) (3) 

is precisely the subgroup of U (4) that conserves the degree of entanglement of 
the states its acts on linearly. That is, 7 G exp(t)) iff V0 G 

H n < 8 ) T-ie where f^ e (0) = 2 \z\ z 4 — z 2 z%\ and 0 and the z 3 are defined in (jT]). 
This Lie group is isomorphic to U( 1) x 577(2) x SU( 2) and comprises all 
matrices of the form C e ® Cn, where £ e , C n are 2 x 2 unitary matrices acting 
separately on the indivdual spin substates. 

Proof. The precise condition for conservation of is readily shown to be 
equivalent to 7 = e H where H, an “inhnitessimal S> e conserver”, fulfills: 

H T Q + QH = i(j>oQ for some 0o £ K; Q = Q 1 — j ® j; 

H = i(j) 0 1 (8) 1 + q 13 i <g> 1 - pi 3 j <g> 1 + k< 8 > 1+ (4) 

912 1 ® i -P12 1 ® j + 1 ® k; ( 3 j, 0o, Pjk, Qjk € K 

and the claimed isomorphism follows straight away from the definition of 
0 . □ 
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2 Lie Theoretic Grounding 


I now turn to the theorem alluded to in the abstract. The full solution of 
several classic applied mathematics problems rests on this proof; such prob¬ 
lems as (i) the Nelson Parking Problem [71ITO] . where the exp (tj Xj) are the 
changes in a car’s configuration (position and orientation on Euclidean M 2 ) 
wrought by driving forwards or backwards in a car with a constant path 
curvature (constant steering setting) and the smallest Lie algebra contain¬ 
ing the Xj is the whole of e(2) or (ii) Montgomery’s description[5] of the 
torque free falling cat as a fibre bundle with the space of cat shapes S as 
the base space, the space of cat’s orientation in an inertial frame (described 
by an element of SO( 3)) as the fibre and the structure (gauge) group is 
some Lie subgroup of SO( 3), depending on what the smallest Lie algebra 
generated by the shape shifting moves the cat can make is. The connexion 
on the bundle arises from computing the shift in the cat’s orientation that 
must accompany a shape shift so as to conserve the cat’s total angular mo¬ 
mentum. One last problem|LIj. less noteworthy than these but worth citing 
here because I personally have formerly dealt with is the proof that a finite 
system of N planar single mode waveguides can realise any transfer matrix 
in U(N ) even though there is only nearest neighbour coupling between the 
waveguides. Mathematically, this means that the Xj span only tridiagonal 
skew-Hermitian N x N matrices, but the smallest Lie algebra containing 
these is the whole of U(N). Therefore (since U(N ) is connected), a finite 
sequence of nearest neighbour coupled waveguides can realise any transfer 
matrix in U(N) by dint of Theorem [0 A second reason for mentioning this 
particular application is that, even though Theorem [D and the methods of 
the present paper give a solution that can be broadened to any Lie group 
and algebra, in the case of nearest neighbour coupled waveguides there is a 
particular solution which, in my opinion, is superior to my own work for some 
applications and that is the algorithm presented in [Qj. Here a sequence of 
transfer matrices of 2 x 2 symmetric coupled waveguides fore-multiply a gen- 
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eral U(N) member in a procedure somewhat like Gaussian elimination. The 
unitarity of the transfer matrix means that rows and columns are reduced to 
one, unit magnitude element at once (instead of needing reduction separately 
as in Gaussian elimination), thus N — 1 basic steps reduces the problem of 
synthesis of a general element of U(N) to the simpler problem of a realisation 
of a general element of U(N — 1), thence, on application of N — 2 basic steps, 
to a member of U(N — 2) thus the algorithm, inductively, will decompose 
the general U(N) to a realisation as a finite sequence of symmetric couplers 
and phase delays. 

I now state and prove in detail the first theorem cited in the abstract; in¬ 
deed this is a slightly more general result, holding for any set of analytic paths 
a : R —* (3 passing through the identity instead of a set of one-parameter 
groups {exp(tX,)| Xj 6 g; i G R}. Results similar to this one seem often to 
be tacitly assumed in the literature but explicit proofs seem to be few. Many 
texts justify something like the following theorem with a glib appeal to the 

Trotter product formulas 

2 

exp([XF]) = Jim (exp (f) exp (£) exp (-f) exp (-f))” and exp(X + 
Y) = lim (exp (—) exp (—))". However, an argument grounded on these 
formulas only establishes a weaker version (stated in [13]) of the result in 
the abstract, to wit, that the closure of the group generated by the e rXj is 
exp(fj), where f) is the smallest Lie algebra containing the Xj, so that a finite 
product of terms of the form e rXj can be arbitrarily near to any element of 
= (exp(f))). This is not a restriction or weakening if the smallest Lie alge¬ 
bra concerned is that of a closed (and topologically embedded) subgroup of 
the whole, but this of course is not always so. In this latter case, one can then 
strengthen the weaker result to be equivalent to the result proven below for 
GL(N, R) through the result that every linear Lie group is Lie-isomorphic (if 
not the same as) to a closed Lie subgroup of GL(N, R) as shown in |2j. For 
example, the irrational slope one parameter subgroup of the 2 -torus group 
U( 1) x t/(l) is not closed, but it is isomorphic to (R, +) and the latter is a 
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linear Lie group and closed as the subgroup of GL( 2, R) comprising upper 
triangular 2x2 matrices with ones along the leading diagonal. The follow¬ 
ing proof avoids all this complexity to get the stronger result directly and 
constructively; thus can guide the general algorithms of §3J 

Theorem 1 (Steering a Lie Group Path Without the full rank Tangent 
Space). Let 0 be a connected Lie group, g its Lie algebra. Kit the appropriate 
sized nucleus A/jd with geodesic co-ordinates so that the co-ordinate map is 
A = log : Afid —> 0 and now let Oj : [—1, 1] — > A 4i; cx.,(0) = id; d r <7j (r) | r=0 = 
Xj\ j — 1, 2, ■ ■ ■ , M be M C u paths (in the geodesic co-ordinates) through 
0 with aj(—r) = (Jj(r)' 1 . The Xj may or may not span the whole of g; the 
typical situation, and main point of this theorem, is where they do not. 

Let furthermore 1} be the smallest Lie algebra containing the Xj ; otherwise 
put: 1) is the intersection of all Lie algebras containing the Xj and is the set 
of all entities that can be gotten from the Xj by finite sequences of linear 
(scaling and addition) and Lie bracket operations. Then for every I 6 [) 
there is a finite number of terms product of the basic paths defined by: 


R 

a : [-1, 1] -)• 0; cr(r) = ]^[ a k (a k r) (5) 

k =1 

such that the tangent to the path a at the identity is d T cr(r )| r=0 = X. 


Corollary 1 . From Theorem [TJ and from HU it follows that every member 
of exp(f)) and thus every element of the (possibly non topoligically embded- 

OO 

ded) Lie subgroup S) — (J exp(f)) fc corresponding to t) C 0 under the Lie 

k =1 


Q 

Correspondence [T 2] can be realised as a finite product of the form cr J ( fc )(r/ c ) 

k =1 


i.e. as a finite product of terms of the form aj{r k ). 


Proof of Theorem\T\ Given two C u path segments ax, cry '■ [—1, 1] —* 0 
through 0 with ax( 0) = cry(0) = id, ax(—T) = (ctx(t)) -1 , ay(—r) = 
(cry(r)) -1 and with tangents X = d T crx(r)| r=0 , Y = d r OY(r)| T=0 to the 


identity, we first show how to realise a path with tangent [. X , Y] to the 
identity as a finite product of these paths. We assume, of course, that [A", Y] 
must be nonzero and linearly independent of X and Y (otherwise there is 
nothing to prove: the smallest Lie algebra containing X and Y is the vector 
space spanned by them). 

Thus let: 


N - \ 

tX+ V /,,(>) .V, 

tY + E Qj(T)Xj] 

where the Pj(r) and qj(r) are analytic functions of r comprising only second 
and higher powers of r and we must allow for the possibility of all the Xj 
that span g being present as terms multiplied by second and higher powers 
of r. So now we consider the tangent at the identity to the family of C u 
paths parameterised by the parameter s as follows: 

a s : [-1, 1] ->• 0; a s (r) = a x (s ) cr y (r) a x (~s ) (7) 

so that the tangent to a family member at the identity as a function of the 
parameter s is, by elementary calculation: 


crx(r) = exp 
cry(r) = exp 


T(s) 


exp ^ad(Af) s + Pj ( s ) Y 

N—l 

(1 + Po(s)) Y + (s + pi(s)) ad(X) Y + pj(s) Zj 

3 = 2 


where we have absorbed all the higher powers s 2 ad(A") 2 Yj 2!, s 3 ad(X) 3 y/3!, • • • 

N -1 

into the sum ^2 Pj{s) Zj with modified analytic co-efficient functions pj(s). 
j =2 

Here the Zj e g are vectors which are mutually linearly independent and also 
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linearly independent from either Y, [X , Y] or any linear combination of these 
last two. We have allowed for the possibility of the vector ad(X) Y showing 
np in the high order powers of s by adding pi(s) to the multplier s of the 
term ad(A) Y; likewise the p 0 (s) term absorbs any Y component from high 
order powers. The co-efficient, functions pj(s), including different powers of 
s, are all different and so are linearly independent on any interval. Therefore, 
we can choose m discrete Sj E [—1, 1] to yield m linearly independent tan¬ 
gent vectors T(sj) where m < N is the total number of linearly independent 
vectors in the set {Y, [X , Y], ad(A") 2 Y, ad(A) 3 Y, • • • }. Therefore, through 
Gaussian elimination, we can linearly combine the tangent vectors T(sj) so 
that their linear combination equals the term [.A, Y]. That is, we can fold 
superpositions weights otj such that: 

m 

[X, Y] = J> fc T( Sfe ) (9) 

k =1 

Therefore, since the Lie group’s identity tangent space is a linear space, 

m 

the path cr(r) = crx(sk) cry(afc r) crx(—Sk) has the tangent [A, Y] to the 

k=\ 

identity. Note that in most cases, not all the vectors spanning the Lie algebra 
g are present in the sum in (]8j), so that the sum in (J9]) contains fewer (m) 
than N , terms. It may even be that none of these terms are there aside from 

Y and [ X , Y], but the point is that these last two are guaranteed to be in ( 151 ) 
with the nonzero weights shown there, so that the procedure above can be 
summarised as: T(s) includes Y, [ X , Y] and some other linear independent 
vectors in the Lie algebra with linearly independent co-efficient functions of s. 
Therefore, there exists a linear combination of the some T(sfc) that sums to 
any of these linearly independent vectors, thus, in particular we can retrieve 

Y or [ X , Y] in this way. 

Given we now have a set of C u paths with tangents Xj, we can build 
paths with a tangent to the identity equal to any linear superpositon of the 
Xj by the method of ma and we can build a path with a tangent to the 
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identity equal to the Lie bracket of any pair of these paths with the procedure 
described above. Therefore a finite sequence of operations comprising the 
above procedure together with linear superposition operations by the method 
of HP can realise a C u path with any tangent inside some neighbourhood H 
of 0 in the smallest Lie algebra () containing the X y 

Having built paths whose tangents span 1), we can, arguing as in m ,realise 
any member of exp('H) on a finite product of these spanning paths. From 
there we can see that a finite product of such paths can realise any member 

CO CO 

of the connected Lie group f) = U exp('H) fc = U exp□ 

k=1 k=1 

Sifting carefully through this proof, one can specialize it to the paths 
a 3 (r) = exp(Xj) and at the same time reword it to make it almost explicitly 
constructive. 

Algorithm 1. Let 0, A, g and Aid be as in Theorem []] and now specialise 
{<jj : [—1, 1] —> Afid] o~j(r) = exp (r Xj)}^ 1 where are linearly inde¬ 

pendent and let t) bethe smallest Lie algebra containing the X 3 \ aside from 
the specialisation to exp, all definitions are wholly analogous to Theorem CD 
Then the following algorithm terminates, and it does so precisely when it has 
found a basis for the whole of f). 
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o := lengthO := M; done := False; 5 = 0; 

while not done do { 

for each j e {1, • • • lengthO} do { 
for each k e {j, ■ ■ ■ , lengthO} do { 

X : = 0j ; Y := Ok] 

U : = maximal linearly independent set from Y, ad(A) Y, ad(A") 2 Y, ■ ■ 
n : = length (W); 

5 := {si, s 2 , ■■■ , Sn} C M 9 e siad(x) y are ]i liear ly independent; 

} 

A := {e s ‘ ad W7| Si G 5}; 

Cnll any member of set A which is linearly dependent with set o; 

Append culled A to set o; 

} 

if length(o) = lengthO then done := True;else lengthO := length(o); 

} 

return o; 

Here the set of vectors being built is the variable o. 

Proof that Algorithm^ terminates and builds the whole of fp The only not- 
explicitly-constructive part of this algorithm is the finding of 
S := {si, s 2 , ■■■ , Sn} C M 9 e s iad(X)y are li near ly independent inside 
the innermost for loop. But since exp(s ad(A) Y) is a linear combination of 
members of: 

N 

U = \J (ad(A) fe Y} = {Y, ad (A) Y, ad(A) 2 Y, ad(A) 3 ¥,■■■} (10) 

fc =0 

the proof of Theorem [2] straight after (J9]) shows that (i) we can find such a 
set by Gaussian elimination and ( ii) that vector space spanned by this set 
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is precisely the same as the vector space spanned by U. So by dint of this 
existence, we can assume that we have some algorithm that implements this 
one step. 

Now witness that the inner for loops add only vectors linearly independent 
from each other and linearly independent from all members of the set o. Since 
there are at most only finitely many, N of these vectors, the algorithm clearly 
terminates. Hence span(o) C f) at all times. 

Now consider the while loop straight before its very last pass. In this 
pass, as just shown above, it calculates sets of vectors that span every set of 

N 

the form [J {ad(X) fc H} for every pair XY E o; in particular is calculates a 
k =o 

vector space containing [ X , Y] for every pair X, Ye o. Since the algorithm 
terminates, this last pass cannot add new vectors to o. Hence t) C span(o) 
when the algorithm terminates. □ 

I have shown the existence of appropriate S at each step within the in¬ 
nermost for loop, so for the purposes of this paper I shall simply assume that 
we have some “black box” that can reckon appropriate Sj for us. Simple 
numerical conjugate gradient optimisation will be used in jj3j The possible 
choices are highly non-unique: almost any randomly chosen different real 
numbers will work. However, the choice of the Sj does bear strongly on “how 
linearly independent” the calculated basis for 1) turns out to be, i. e. it bears 
on the condition number of the array of linearly independent vectors. This 
condition number in turn will bear on how many basic magnetic field pulse 
sequences will be needed in §3] and with what pulse times. Poor condition 
numbers mean long and complicated pulse sequences comprising long linger¬ 
ing pulses, and therefore this means that in practice the pulse sequence is 
likely to be highly sensitive to errors in the experimentally applied magnetic 
fields or pulse times. In some cases, one may be able to use the Cartan 
decomposition to find a generalised “polar” decomposition for SU( 4) thus 
achieving a better or minimal sequence of evolution operations as done in 
[3], but this last method’s applicability depends on the Cartan pair’s Lie 
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subalgebra 1 C g and its orthogonal (with respect to the Killing form on g) 
complement £-*- being directly available, as operators from exp(l) and exp(t ± ) 
to the experimenter as time evolution operators; in general, however, at least 
one of the available vectors of the form found in Algorithm [Tj will “straddle” 
the Cartan pair (have nonzero components in both t and fA). The finding 
and optimisation of a systematic algorithm, grounded on the Cartan decom¬ 
position or otherwise, for calculating S is therefore a high priority for future 
work on the procedure of (J3J 

Once we have constructed a full basis {Z .,}) 1 ™ 11 for f) with Algorithm 
[U and given any 7 £ S) = (exp(f))) within the connected Lie group A = 

OO 

(exp(f))) = (J exp(f)) fc can then, by Corollary Q] find real numbers 77 & so 
k =1 

jo / dim h / , \ \ 

that 7 = n n exp IrjkZk) ) for some finite jo- There are two more 

j =1 V k=i ' ') 

steps to this construction. Firstly, we need to decompose 7 as 7 = Id'/Li e ' ij 
for Yj £ f). If the Lie group A is compact, for example, then it is the 
exponential of its own Lie algebra we can always find a (nonunique) Y £ t) 
such that 7 = e y . But a general, connected Lie group ( e.g. SL( 2, C) or 
SL( 2, M)) is strictly bigger than simply the exponential of its Lie algebra: 
there are some elements which are not exponentials of the algebra. In )J3] 
we deal with the compact group S'L/(4), so this first step is not needed and 
will not be further considered in this paper, but in general a decomposition 
algorithm must be found. For example, in a connected, semisimple Lie group, 
the Iwasawa decomposition |Tj will achieve this goal. 

Thus I now assume that we must realise 7 = e 5 where Y £ t). Since one 
can write 7 = e 1 = exp (7-)" for any n £ N, there is an n big enough that 
exp (h) g /C, where /C is any arbitrarily small nucleus. We choose the latter’s 
size right and then the Wei-Norman equations^4] can be used to calculate 
the canonical co-ordinates as follows: 

Lemma 2 (Wei-Norman). Let 0 be an N dimensional connected Lie group, 
g its Lie algebra, {Aj-}A =1 C g vectors spanning g and /C C 0 a nucleus 
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small enough to be labelled by the canonical co-ordinates of the second kind, 


N 


i.e. K, C < n e TjXj \ Tj E [—1, 1] >. Let 7 = e x E /C, then the canonical co- 
3 = 1 


ordinates {rj} x =1 of the second kind for 7 = e x where X = x\ X\ +x 2 X 2 +- ■ ■ 
are calculated by solving the TV-dimensional nonlinear differential equation: 


d t T = M(ti, r 2 , • 


M(ri, t 2 , •••) = 


( X\ ^ 

x 2 

X3 

\ / 


/1 

^ ^ e nad(ii)j ^gTi ad(Vi) gT2 ad(V2) 


\ 


^ e n ad(Xi) 


. gUv-l ad(Xjv_i) 


TV 


V ; 

subject to the initial conditions r = 0 at t — 0 . 
co-ordinates are the values in the vector r at t — 1 , 


in) 

The required canonical 




TV 


i.e. X = Yl e' J 
3 = 1 


i(l) 


Proof. See j!4j for details. The essential trick is to break the derivative 
d t e n h) Xl e T2(t)x 2 ... U p j n t 0 summands by the product rule, and then to 

g1~m (d V m ^ ... g UV (d -V/V a g 

gn ad(Vi) gT 2 ad(A _ 2 )... g r 2 ad(v m _i) _d_ r (t) X e T1 h) e T2 ^ X2 • • • e T2 ^ Xn bv the 
Wei-Norman “shuffle” trick. One then uses this to write the differential equa¬ 
tion describing the path [— 1 , 1 ] —> /C; t e tx in the form 

X e tx = Xj + f 2 e Tl ad ^i) X 2 + f 2 e n ad (*i) e T2 ad (^ 2 ) jy 3 • • • j e tx , cancels 
the terms and inverts the matrix to get (ITT]) . The matrix to be inverted 
in (ITT]) is such that its j th column is the j th column of the square matrix 

gTi ad(A'i) gT 2 ad(A 2 ) . . . gTj-i ad(V,_i) ^ w } len the g-r fc ad(X7) are wr if t en aS matrices 


rewrite each term of the form e T1 ^ Xl e T 2 h ) A ' 2 • • • -L 
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with respect to the basis {}^ =1 ). The matrix to be inverted is the iden¬ 
tity matrix at r = 0 and its determinant is a continuous (indeed analytic) 
function of the Tj, therefore there is a nonzero interval — e < t < +e for 
e > 0 wherein the matrix is invertible and so (IIip is a Cauchy initial value 
problem fulfilling the conditions of the Picard-Lindelof theorem, i.e. that 
the vector function on the right of (fill) is Lipschitz continuous function of 
the Tj. Therefore, as long as X G g is near enough to 0, the unique, existing 
solution of the Cauchy initial value problem defined by (1X11) must define the 
co-ordinates for X at t = 1. □ 

and this lemma yields the algorithm to realise any 7 G 0 : 

Algorithm 2. Let 0, N and g be as in Theorem CD g be spanned by in¬ 
dependent vectors {Xj}^ =1 and suppose further that 0 = exp(g), i.e. that 
every group member can be written as an exponential of a Lie algebra mem¬ 
ber. Let 7 = e x ; X 6 g be an arbitrary member of 0. Then the following 
algorithm: 
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n := 1; det := 1; t := 0 ; r = 0; 

Y := X- 

for t e [0, 1] do { 

Calculate r(f) by Wei-Norman equations for Y ; 
det = det(W _1 ); 

if det < e break; 

} 

if det > ethen return t, 1; else 
{ 

n : = Floor (|) ; 

Y := 

} 

for t e [0, 1] do { 

Calculate r(t) by Wei-Norman equations for Y; 

} 

return r, n 

terminates and outputs the canonical co-ordinates for = exp (^-) in r and 
the power n that one must raise 7 n to to realise the original 7 = e x . 

Proof of efficacy and convergence of Algorithm. Clear from Lemma [ 2 j □ 

The for loop is shown over a continuous iterator variable t, so that, as 
written, the loop must run Mi times! In practice, the loop is a discrete 
numerical solving of the differential equation and it breaks, unsuccessfully, if 
the determinant of the matrix drops below a critical threshold value e. Since 
the Wei Norman equations will succeed (barring some numerical difficulty) 
in finding co-ordinates for near enough to the identity, either the first for 
loop (with n = 1 ) will return the correct co-ordinates or the second for loop 
will do so for A with n now calculated to be high enough to bring exp (^) 
near enough to the identity to belong to the nucleus /C of the statement of 
Lemma |2j 
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This ends our Lie theoretical background discussion and I now turn to 
the synthesis of time evolution operators in SU(4) for coupled two qubit 
quantum systems. 

3 Coupled Spin State Manipulation 

To further analyse the control that the experimenter wields in applying se¬ 
quences of magnetic held pulses to the coupled spin system of (JU I define 
the following 4x4 traceless, skew-Hermitian matrices: 


X = X 0 + kK Y — Yq + kK Z = Z 0 + kK 
k x = | (-j <£> k + k <g) j) k Y — I (k 0 i - i 0 k) 
k z = 5 (-i 0 j + j ® i) k xx = -f(j®j + k 0 k) 

K Y y — ~\ (i 0 i + k<2) k) K zz = (i 0 i + j 0 j) 

L x = (i 0 1 + 1 0 i) Ly=ja®l + l®j) 

L z = -|(k ® 1 + 1 ® k) (12) 

k X Y = — 2 ( 7n + 7e ) (7ni® j +7ej ® i) - f (i® 1 - 1 ® i) 
k Y z = - 2(7 ; +7e) (7nj ® k + 7ek ® j) - f (-j ® 1 + 1 ® j) 

kzz = 2(7 ; +7e) (in k ® i + 7e i ® k) - f (k ® 1 - 1 ® k) 

X KK = |(i®l-l0i) Y kk = \ (—j 0 1 + 1 0 j) 

Z K k = |(k0l — 10 k) 

Lemma 3. Given the dehnitions of (J 2 ]) and dT 2 |) and the fixed, strictly pos¬ 
itive real number k we have: 

(a) The smallest Lie algebra containing {B x X + B y Y + B z Z ± K\ Bj E M} 
is the whole of su(4) whenever at least one of the positive constants y n , y e 
are nonzero; 

(b) Suppose that at least one of the positive constants 7 n , y e are nonzero and 
that 7 n ^ 7 e . Let X (6 , 0, a) = sin 6 cos 0 X + sin 6 sin (j) Y + cos 6 Z + a K. 

def 

Then the smallest Lie algebra containing b(d, 0) = {X(9, 0, a)| a E M} (for 
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fixed 9, 0 e M) is the five dimensional Lie algebra 
b(0, 0) = span R ({X e ^, K, ad(Xg j<j> ) K, ad(Xg !(f) ) 2 K, ad(K ) 2 X e ^}) (where 
we write Xg^ = X(9, 0, 0)) and exp(b(0, 0)) is a Lie group isomorphic either 
to SU ( 2 ) x (U( 1) x R) or SU(2) x (f/(l) x 17(1)) = S'L/( 2 ) x T 2 . 

(c) If 7n = 7e 7 ^ 0, then the Lie algebra b(0, 0) in statement (b) shrinks to 
the four dimensional b (9, 0) = span R ({X< 9 ;( y Ji, ad(Xg j( p) K, ad{Xg^) 2 K}) 
and exp(b(0, 0)) is a Lie group isomorphic to 517(2) x 17(1); 

(d) If at least one of 7 n , 7 e are nonzero and 7 ^ 7 e , then the smallest Lie 
algebra containing both b( 0 i, 0 i) and b( 0 2 , 02 ) for at least one of 0 i 7 ^ 0 2 , 
9 1 7 ^ 0 2 is the whole of su(4) and its exponential is the whole of 517(4); 

(e) If 7n = 7e 7 ^ 0 then the smallest Lie algebra containing both b(0 1; <f>i) 
and b(0 2 , 0 2 ) for at least one of 0 1 7 ^ 0 2 is no longer the whole of su(4) 
but is instead a nine dimensional Lie algebra with su(3) as an ideal and its 
exponential is 517(3) x 17(1) (naturally with 517(3) a normal subgroup). 

Proof, (a) The fifteen matrices X, Y . Z. K x ■ Ky ■ K Z , K XXl X YY , K Z z , Lx, 
L y . L z , K X y 1 K yz , K Z x are readily shown to be unconditionally linearly 
independent. Indeed if they are vectorised into their 16 real-component su¬ 
perposition weights over the basis {ly ® u*| j, k e 0 • • • 3} and u 0 = i 1, ui = 
i, u 2 = j, u 3 = k and built, together with the same vectorisation of i times 
the 4x4 identity, into a square matrix of 16, 16-component row-vectors, then 
the determinant of this array is always — ^ as long as at least one of y n , y e 
are nonzero. These fifteen are derived from the algebra members X, Y, Z, K 
by the bracket relationships: 
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K x = ± 
Kv = \ 

k z = \ 

Kxy = l 
Lz — | 


K, X 
K, Y 

k, z 


Xrk — 


Yrr — 


Zrk — 


1 

2 

k x , Y 

Kyz = | 

1-1 

0 

7- 

i-1 

Lx — \ 


K, K x 

k, ky 

k, k z 
k Yl z 
k Y , kz 


K X x = h 


Kyy = 


k zz — 


1 


K x , X 
Ky, Y 

k z , z 


+ « k KK 
+ « ykk 
+ K Zrk 


Kzx = i 


^ZX 


Ly = \ 


Kz , X 


K z , K 


x 


(13) 

Thus (a) is proven, but we note in passing that an alternative linearly inde¬ 
pendent set is as above but with L x , Ly, L z replaced by X k K, Y kk , Z kk . 
However, in this case, the determinant stated is — 4 ( 7e ~ 7n 1 and vanishes 
when 7 n = 7 e . So, although we use this latter set for our realisation below 
(it is simpler), it cannot be used for coupled spin systems when y e ~ 7 n . 

(b) Without loss of generality, one can align the co-ordinate axes to the 
field to show by simple computation that e.g. {Z, K, K z , k Z z, Z KK } is the 
claimed smallest algebra, that {Kz, Kzz, Zrr} is an ideal of this algebra 
and that the ideal is isomorphic to su(2). For the three separate ideals 
spanned by the K x • • •, k Y ■ ■ ■ and K z ■ ■ ■ matrices, it is readily shown 
that there are two seperate central elements for the exponentials of these 
ideals, so that the fundamental groups of these exponentials are the same as 
those of 577(2), not 50(3). Lastly, [Z, K] lies inside the ideal, so therefore 
the corresponding cosets commute in in the factor algebra and the quotient 
group must be a commutative Lie group. e Kt is, as we have seen, periodic, 
thus the one parameter group defined by K is 0(1). The one parameter 
group defined by Z is 0(1) if 7 e / 7 n or 7 n / 7 e is a rational number, otherwise 
it is the noncompact group (M, +), analogous to the irrational slope one 
parameter subgroup on a 2-torus. 

(c) Readily proven as for (b). In this case, it is {K z , K Z Zi Z} is the ideal 
exponentiating to 50(2). 

(d) Without loss of generalness, we can assume X, Y as the two directions: 
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if there are two linearly independent directions, then these span a plane so 
that we can choose orthogonal vectors spanning the plane. Then we derive 
Kxy, K\, Ky from I\. X, Y as in (fl3l) and then form: 


[Kxy, X} + k(^ 


V ln+le 


[K x , Ky} + 


(7n~7e) 

7n+7e 


K + 2 k 2 K x + 2 


7e7» 

(7e+7n) 2 


K 


Y 


_ K(7n~7e) 2 , 

'In+le 

(14) 

so that, as long as y e ^ y n , a finite sequence of linear and bracket operations 
yields the vector Z, whence the generation of the whole algebra proceeds 
exactly as in (a). 

(e) As can be checked by straightforward but tedious calculation, the smallest 
algebra generated by A', X, Y when y e = y n is 
span M ({i7, X, Y, K x , K Y ,K X x, Kyy , K X y, The following linear 

combinations: 


At = a; 




K 


YY 


X 2 = 2(k XY -kX + k 2 K) A 3 = 2 Li 


\ _ X-kK+Ky 

Xi ~ -75- 

x _ X-nk-ky 
Ae “ 72 


\ _ Y-kK-I<x 

As ~-75- 


(15) 


A 7 — 


Y-kK+Kx 

72 


\ _ Kxx+Kyy 

As - —73— 


span an eight-dimensional ideal, thus exponentiating to a normal subgroup 
and the factor group is the one parameter group e Kt , which is periodic, thus 
isomorphic to U( 1). The eight vectors in (1T51) have the same structure con¬ 
stants as those of the Gell-Mann matrices (divided by i, since the physicist’s 
convention is work with the algebra of Hermitian matrices and to form ex¬ 
ponentials as e lH ) and thus of su(3). The basis in (1131) begets the diagonal 
Killing form Tr(ad(Aj) ad(Afc)) = —12 5j k with A x , A 2 , A 3 spanning a copy of 
su( 2 ), As commuting with any of these three, Ai, A 2 , A 3 , As spanning the Lie 
subalgebra in a Cartan pair of a Cartan decomposition and the other four ma¬ 
trices spanning the orthogonal complement in the Cartan pair. Now we must 
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decide whether the Lie algebra spanned by (fT5]l exponentiates to the whole 
of SU(3) or a discrete projection thereof (PSU(3)). For SU( 3), the group’s 
three-element center {exp 1) l| k G Z} is contained within the one pa¬ 
rameter subgroup corresponding to j exp ^(3 A 3 + 2 \/3 A 3 ) r j re R j . For 
the matrices in (1T5]) . it is also readily shown that there are three distinct, 
central elements of the form exp ^(3 A 3 + 2 a/3 A 3 ) I/p j ; k e Z, so the fun¬ 
damental group and Lie algebra of the exponential are the same as that of 
577(3) and the exponential is therefore isomorphic to SU(3), not a projection 
(PSU( 3)) thereof. □ 


Practical consequences for experimenters The practical outcomes of 
Lemma [3l given Theorem Q] and Algorithm [TJ are (i) that an experimental 
kit which can impart variable magnetic fields of either sign in three, lin¬ 
early independent directions to a coupled spin system can always impart any 
time evolution operator from SU (4) to the two coupled quantum spin state, 
(ii) that a simplified experimental kit imparting magnetic fields in only two 
linearly independent directions can do the same (at the expense of more 
complicated magnetic field pulse sequences) as long as q n / 7 e , so that for 
coupled spin systems with similar or the same gyromagetic ratios for the two 
coupled particles, all three independent magnetic fields will be needed for 
good controllability of the time evolution operator. 

Lastly, 1 now give an example of the algorithms’ application for the cou¬ 
pled 31 P nucleus-electron system. We first need to calculate a basis for 
su(4) with Algorithm [H To this end, let us define our “unit” magnetic field 
strength as B u = lOmT and our “unit” time interval as r u = 100ns. Linder 
these unit conditions, we define K v = 5.8765 A", Xu = B 0 (y n + 7 / t u X 0 = 
27.98723 X 0 , Yu = 27.98723 Y 0 and Z v = 27.98723 Z 0 . By imparting a 
magnetic field (possibly none) B Xl B y , B z e R units for a time interval r 
units, we can now directly impart any time evolution operator of the form 
exp ((B x , Xu + By , Yu + B z Zu) + Ku)r) for r > 0. As discussed in we can 
achieve a term of the form —Kut for t > 0 as +Ku ( nr p — t ) where n e N 
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and T p ~ 2.14 in our 100ns time units is the period of e Kut . 

By examination of the steps of TheorernQ] and Algorithm [TJ we can im¬ 
part time evolution operators belonging to the one parameter subgroups 
{e rHj \ t G M} corresponding to the following Lie algebra members Hf 


H\ = (h Xu + Ku) n 

H 2 = ( b 2 Yu + Ku)T 2 

H.i = (63 Zu + Ku) r 3 

H 4 = t 4 Ru 

H 5 = Ad(exp((&5 Xu + Ku) r 5 )) H 4 

H 7 = Ad (exp (AT Ty)) Hi 

H 6 = Ad(exp((6 6 X v + K v ) r 6 )) H 4 

H 8 = Ad(exp((6 8 % + Ku) r 8 )) A 4 

H w = Ad(exp(Ku t 10 )) H 2 

Hg = Ad(exp((6 9 % + K v ) r 9 )) H 4 

Hu = Ad(exp((6 n X v + K v ) m)) H 4 

H \2 = Ad(exp((6 12 Xu + Ku) t 12 )) H 4 

H13 = Ad(exp(( 6 i 3 X v + Ku) n 3 )) Ad(exp((ci 3 Y v + I<u) ^13)) H 4 


Hu = Ad(exp((&i 4 Yu + Kjj) tu)) Ad(exp((ci 4 Z v + Kjj) <p 4 )) H 4 
H lb = Ad(exp((6i 5 Y v + K v ) n 5 )) Ad(exp((ci 5 Z v + Kjj ) <Ji 5 )) H 4 

(16) 


We can here use the second basis (with Xkk, Ykk, Zkk ) discussed in 
the proof of Lemm [3j part (a). We can understand why the forms in (fl6l) 
do indeed yield an independent basis by looking carefully at the proof of 
Theorem [2] and Algorithm |T| Firstly, //].•••, 7/ t are linearly independent. 
The forms of H 5} • • • , iL 6 , H 8} H 9 , H u , H V2 are the forms in (JTO and (f 8 j) 
with a x set to paths of the form exp(s Hj) for j G 1 • • - 3 and Y set to ay 
set to exp (t K u ) with different values of .s in (J8J) . Since we have seen that 
Xu and Ku generates a Lie algebra of dimension 5, there can only be two 
linearly independent values of the form of H 5 and Hq and one of the form H 7 
realised with different values of sin (J8j). Likewise for H$, Hg, H\g. We realise 
two further linearly independent values Hu, H \ 2 likewise but when try to 
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get a further independent value of the form Ad (exp (A^ t 7 )) H 3 we fail. This 
is because we have already included Ku in our list of linearly independent 
vectors, and then Zkk is linearly depenent on all the others through the 
relationship X KK + Y KK + Z KK = —2 AT. Lastly, the forms of iL 13 , H Ul i7 15 
give us the last three linearly independent vectors because we have seen in 
Lemma (lT 6 j) that ad(X) ad(F) A",ad(F) ad(Z) K and ad(Z) ad(X) A' yield 
the last three members of the su(4) basis. 

Now we must choose the numerical values of the constants in (fl 6 |) . This 
was done numerically by vectorising all the fifteen Hj (calculating their com¬ 
ponents with respect to the basis {in,- <gi u k \ j, k e 0 • • • 3} with u 0 = i 1, u 7 = 
i, u 2 = j, u 3 = k), assembling these components into a 15 x 16 real matrix 
and then calculating the singular values of this matrix. The condition number 
(ratio of maximum to minimum singular value) measures it how linearly in¬ 
dependent the basis is and therefore how well it is likely to work in Algorithm 
[2l That is, poor (large) condition numbers mean higher values of n output by 
Algorithm [2] and therefore much longer preparation magnetic pulse sequences. 
So the constants of (1T6l) were optimised with a simple hill climbing procedure 
(repeatedly working through the list, optimising each constant individually) 
to find the constants yielding the minimum condition number for the basis. 
The following constants yielded a maximum singular value of 15.82 and a 
minimum singular value of 1.70, i.e. a condition number of 9.3: b\ = 2.003, 
Ti = 0.151, b 2 = 1.5155, r 2 = 0.176, b 3 = 1.958, r 3 = 0.118, r 4 = 1.109, 
b 5 = 0.3015, r 5 = 1.021, b 6 = 0.5195, r 6 = 0.910, r 7 = 0.215, b 8 = 0.1925, 
r 8 = 0.931, b 9 = 0.222, r 9 = 0.926, no = 0.2005, b u = -0.167, nl = 0.9825, 
&i 2 = 0.394, T 12 = 0.9255, b Vi = 0.198, n 3 = 1-017, ci 3 = 0.178, <q 3 = 0.9855, 
&i 4 = 0.344, ri 4 = 1.000, C 14 = 0.166, <ji 4 = 0.9845, 615 = 0.257,ns = 0.900, 
C 15 = 0.190, ns = 1-377. 

Let us now use this basis to realise the unitary SU (4) matrix 7 = j < 8 > i. 
This has a physical meaning of the change that we should see if we could 
impart a magnetic field in the y direction on the electron, a magnetic field of 


24 


the same magnitude but in the x direction on the nucleus and allow both of 
these spin states to precess, uncoupled and independently, a one-fourth cycle 
about their respective magnetic fields. Clearly there is no simple, obvious 
way of doing this with the fields available to the experimenter. The matrix 
logarithm of 7 is log 7 = —z ^ (i®j + k®k) and its approximate components 
in the basis of (fl6l) are x\ = —0.287893, x 2 = —0.41226, x 3 = 0.178931, 
x 4 = -0.846392, x 5 = 0.248348, x 6 = 0.215918, x 7 = 0.163212, a: 8 = 0.77681, 
xg = 0.143116, X 10 = 0.204211, xu = 0.19128, x 42 = 0.219224, x\ 3 = 
0.517963, x 14 = -0.363032, x 15 = -0.269563. 

We now use a Runge-Kutta discretised approximation of Algorithm [21 
One thousand steps (At = 0.001) is found to be extremely accurate. When 
we try to find the canonical co-ordinates of 7 directly by the first for loop, 
the determinant of the Wei-Norman matrix vanishes at about step number 
166, i.e. at about one sixth the way through the loop, signalling that roughly 
exp (7 log 7 ) is the furthest matrix from the identity that the Wei Norman 
equations can find second kind canonical co-ordinates for. Therefore, the 
second for loop in Algorithm [2l we set n = 7 and find the canonical co¬ 
ordinates of exp (y log 7 ). Now the Wei-Norman equations can be solved 
and yield the following canonical co-ordinates of the second kind: 
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(M£ i 


{ -0.0724497 -0.0480297 0.0459191 -0.0680989 0.0333669 

0.0469188 0.0328527 -0.125384 0.203289 0.0477131 

0.019957 0.0193781 0.0793817 -0.0672283 -0.0156309 } 


15 

exp (i log 7) = 


3 = 1 


g hj Hj ~ /-y 1 ' ^ 


( Til 
0 
0 

V -714 


0 0 714 \ 

722 - 7*4 0 

7*4 722 0 


0 


0 


711 / 


711 = 0.950484 + 0.216942i; 714 = -0.216942 - 0.0495156i 
722 = 0.950484 - 0.2169422 

( 17 ) 

and the rms error of the matix elements y // ^tr((7^_ — 7)! (7^ — 7) is found 
to be 6 x 10 12 , where 71 is the numerical value defined by the numerical 
second kind canonical co-ordinates in (fT7j) and 71 is the numerical result 

7 

of Algorithm [21 i.e. the numerical realisation of our design goal 7 = j <S> i. 
The lengths of the magnetic pulses found here are of the order of tens to 
hundreds of nanoseconds, and the speed of pulse manipulation is thus well 
below the 1GHz upper frequency limits of the contemporary magnetic held 
control speeds cited in [5]. The peak magnetic fields strengths used here 
are of the order of 2mT (e.g. the magnetic held in h\ {b\ X rj + Ku), which 
is 2.003 x 0.287893 ~ 0.58, corresponds to 5.8mT, and the constants h 3 
in the canonical co-ordinates are between 0.1 and 0.3); these are big but 
are quite in keeping with the peak manipulation helds available in today’s 
technology. Lastly, one Wei-Norman cycle, i.e. preparation sequence of the 
form e hlHl e h2H2 ■ ■■e hlsHl5 takes forty six magnetic pulsing/ idling stages. 
The whole preparation therefore takes 322 pulsing / idling stages, each of the 
order of ten nanoseconds. Assuming this example to be typical, the whole 
preparation sequence takes of the order of 5 microseconds, well smaller than 
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the two second coherence time cited in [5], 


4 Conclusion 

In a general Lie group © setting, we have discussed the set generated by Lie 
group members of the form e Tj Xj where Xj e g do not span a Lie algebra and 
given a semiconstructive proof of the fact that they generate, with products 
of a finite number of terms, the (possible non-topologically embedded) Lie 
subgroup generated by exp(f)), where I) is the smallest Lie algebra containing 
the Xj. This proof becomes constructive when the Lie algebra exponentiates 
to the whole Lie group, as happens in a compact group, for example, and 
we have given algorithms and proofs of their convergence in such a case. 
In the case of two qubit coupled spin | quantum states, these proofs show 
that the ability to pulse three linearly independent magnetic fields at any 
field strength in an interval around nought for any positive time allows the 
experimenter to impart any unitary transformation on the quantum state 
space in a finite sequence of these operations. Indeed, when the gyromagnetic 
ratios for the coupled spins are strongly unalike, two linearly independent 
magnetic field directions will suffice. A numerical version of the algorithms 
was run for an example transformation and found to be highly accurate. 
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